Interface Spin-Orbit Coupling in a Non-centrosymmetric Thin-Film Superconductor 
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■ We present a detailed study of the effects of interface spin-orbit coupling (ISOC) on the critical 

field behavior of non-centrosymmetric (NCS), ultra-thin superconducting Be/Au bilayers. Parallel 
field measurements were made in bilayers with Be thicknesses in the range of d = 2 — 30 nm and 
Au coverages of 0.5 nm. Though the Au had no significant effect on the superconducting gap, it 
produced profound changes in the spin states of the system. In particular, the parallel critical field 
exceeded the Clogston limit by an order of magnitude in the thinnest films studied. In addition, the 
parallel critical field unexpectedly scaled as H c m/A oc 1/d suggesting that the spin-orbit coupling 
energy was proportional to A /d 2 . Tilted field measurements showed that contrary to recent theory, 
the ISOC induces a large in-plane superconducting susceptibility but only a very small transverse 
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One of the most fundamental characteristics of a superconductor is the symmetry of its condensate wavefunction. 
Indeed, this has proven to be a central issue in the description of a number of non-conventional superconductors 
such as high-T c systems Q and the ruthenates @> Q. In conventional BCS superconductivity the condensate is 
time reversal invariant and is formed from Cooper pairs consisting of electrons of opposite momentum and opposite 
spin In non-conventional superconductors, however, this simple symmetry can be modified by the underlying 
£N| , crystal structure and/or the symmetry of the pairing interaction. A compelling example of the former is the recently 
discovered heavy fermion superconductor CePt 3 Si, whose crystal structure lacks inversion symmetry CePtsSi 
exhibits a line node gap structure 0,0, 111 which is believed to be, in part, a consequence of strong spin-orbit (SO) 
coupling in a non-centrosymmetric (NCS) crystal symmetry d, 0, llfl | . This has stimulated renewed interest in the 
, possibility of realizing non-conventional pairing states from the convolution of broken inversion symmetry and spin 
orbit coupling, neither of which violate time reversal invariance. Recent theoretical work on the problem has suggested 
that a bilayer configuration, which, by definition, lacks inversion symmetry, consisting of a low atomic mass thin- 
film superconductor coated with a high atomic mass non-superconducting element such as Au or Pt would be a 
model system for studying the effects of SO coupling in a NCS superconductor In the present Letter, 

we present a study of interface SO coupling in thin Be films coated with 0.5 nm of Au via critical magnetic field 
measurements. Contrary to recent theoretical predictions, we find that the SO induced in-plane superconducting 
spin susceptibility is significantly larger than the corresponding perpendicular susceptibility. Furthermore, the SO 
coupling cannot be described in terms of Abrikosov-Gorkov impurity formalism Q in that it appears to be a function 
of the superconducting gap. 

In the experiments described below we use critical field measurements to determine the SO coupling strength in 
Be/Au films of varying Be thickness. The Maki equation [Til Il5| is a useful tool for extracting the spin response of 
the superconductor from the orbital response, particularly when the film is not in the thin-film limit and/or the field 
is not parallel to the film surface. In general the critical field of a thin film is a function of the superconducting gap 
A Q , the film thickness d, the electron diffusivity D, and the spin orbit coupling parameter b. The critical field H c is 
determined by the implicit function |l5| : 
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T c is the critical temperature, and ip is the digamma function, e is a function of the angle between the plane of the 
film and the magnetic field, 

e(0) = L>[2e J ff c sin(6») + \{deH c cos{9) f/h] 

where 9 — corresponds to field parallel to the film plane. The parallel field solutions of Eq. in the d — > limit, 
where orbital effects are negligible, are of particular relevance , 

2 ^ gc " w V2 6/A «l (2) 



= VWA^ > f, (3) 

where Eq. @ is the familiar Clogston spin-paramagnetic critical field 

13 13 EH Note that from Eq. © the critical 
field can be arbitrarily high with increasing b. In contrast, if the Zeeman coupling is neglected then, at any finite 
thickness d < £, the parallel critical field is limited by the orbital term, 



2hbH c1{ _ / 3ft 3 1 

A y m 2 DA d { ' 

Numerous studies of the spin-paramagnetic transition in ultra-thin Al and Be films have shown that these two 
light elements have a very low intrinsic SO scattering rate 0, |3 EH an d are true spin-singlet superconductors. 
Co nseq uently, they make ideal candidates for systematic studies of the effects of ISOC induced by impurity coatings 
fill l22| . Not only does this bilayer configuration break inversion symmetry, but one can treat the ISOC as a boundary 
condition on the spin component of the superconducting wave function. Early measurements of the critical field 
behavior of Al films coated with less than 0.2 nm of Pt revealed that the Pt-induced ISOC can produce significant 
enhancements of the parallel critical field with virtually no effect on the transition temperature [19f. Though the 
temperature and angular dependencies of the Al/Pt critical fields were in reasonable agreement with the Maki formula, 
the overall magnitude of H c u was typically ~ 10 — 30% higher than could be accounted for by reasonable estimates 
of b. 

Recent analysis of the spin states of two-dimensional NCS superconductors, and, in particular, superconducting- 
normal metal bilayers, predicts that ISOC will introduce an anisotropic spin triplet component into superconducting 
ground state 0,[p|- The triplet component is most directly manifest in the spin susceptibility of the superconductor 
X s - Of course, conventional spin-singlet BSC superconductors have \ s ~ at T = 0, but ISOC is expected to lift 
the spin degeneracy, and in the strong coupling 1: X\\ ~ X™/2, xl ~ X™: where xjj is the in-plane superconducting 
susceptibility and \ n 1S the normal state susceptibility. The Clogston critical field given by Eq. is derived by 
equating the magnetic free energy difference between normal and superconducting states with the condensation 
energy of the superconducting state. If \ s 7^ then Eq. @ can be generalized, 

- V2A /y/l-x a /x n - (5) 

This would imply that the parallel critical field of a bilayer, such as the Al/Pt samples of Ref. flij. would never be larger 
than H c u = A //j,b- However, the critical field enhancements observed in those early experiments were significantly 
greater than this upper limit. Since the late 1970's very little work has been done on SO coupling in thin film 
superconductors, though several fundamental questions remain unanswered. For instance, the Al/Pt measurements 
were made at a single Al thickness and over a very narrow range of gap values. Clearly, one would like to determine 
the depth to which the ISOC penetrates into the superconducting film and how the spin component of critical field 
behavior scales with film thickness and gap magnitude. Another important issue is whether or not the anisotropy in 
X s is observable in the angular dependence of the critical field. 

Be/Au bilayer films of varying Be thickness were prepared by e-beam evaporation in a initial vacuum of ~ 0.2/z Torr. 
All of the depositions were made on fire polished glass substrates held at 84 K. First a layer of Be film with thickness 
in the range 2. — 30.0 nm was deposited at a rate of 0.14 nm/s, then a 0.5 nm Au film was deposited at 0.01 nm/s 
without breaking the vacuum. Both the Be and Bc/Au films were found to be very smooth and homogenous, with 
no evidence of islanding or granularity, see Fig. ^ The films were subsequently trimmed in order to eliminate edge 
effects. Resistive measurements were made in a dilution refrigerator with a base temperature 50 mK by a standard 
four-probe lock-in method. The films were aligned with the magnetic field via an in situ mechanical rotator. In the 
data presented below the transition temperatures were defined by the temperature at which the resistance fell to 10% 
of its normal state value and the critical field was determined by the midpoint of the resistive transition. 
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The transition temperature of the homogeneously disordered Be/Au films used in this study are plotted as a function 
of Be thickness in Fig. [5] Films with d < 2 nm are known to display a non-perturbative zero bias anomaly in their 
tunneling density of states, which is associated with the emergence of the Coulomb gap |2^|. As can be seen in Fig. 
121 this is also the critical thickness below which the zero temperature superconducting phase is lost and the electron 
diffusivity goes to zero. In order to make use of Eq. (JIJ it was necessary to fit the thickness dependence of T c and 
D with an empirical functional form. In particular, the solid lines in the main panel and the inset of Fig. [21 are the 
functions: 

T c {d) = T co tanh[(d- 1.35)/1.29)] (6) 
D{d) = L> tanh[d 2 /23.3], (7) 

where T co = 0.68 K and D a ~ 3H/m. D a was determined from films with d > 10 nm using the relation 1/R = 
2e 2 voD d, where R is the sheet resistance and Vq is the density of states per spin of Be. The BCS coherence length 
for a Be film with T c ~ 0.7 K is £ Q ~ 4 fim. For the range of diffusivities plotted in the inset of Fig. [2]the mean free 
path is always l a < 1 nm and the corresponding Pippard coherence lengths are in the range £ = 0.85\/£ o Z o « 20 — 40 
nm. Consequently, all of the data discussed below is in the 'dirty' limit where l -C d < £. 

We have measured the parallel critical field H c \\ of Be films of varying thickness (d m 2 — 30 nm) with and without 
0.5 nm Au overlayers. The Au coatings did not significantly affect T c but did increase the normal state resistances 
by 10 — 50%. Since the transition temperature of the films varied greatly over this thickness range, we normalized the 
critical fields by the superconducting gap. In Fig.[3|we plot 2^bH c u/ A q as a function of the inverse Be thickness at 
60 mK, where we used the relation A Q = 2.1/cbT c , as determined by tunneling density of states measurements |20| . 
The triangular symbols correspond to pure Be films. Note that, except for the thickest Be sample, the critical field is 
independent of d and precisely that of Eq. J2J) , indicating that orbital effects were negligible in films of thickness less 
than ~ 10 nm. In contrast, the normalized critical field of the Be/Au bilayers (circular symbols) is not only a strong 
function of thickness, but it exceeds the Clogston limit by more than a factor of 8 in the thinnest film! Furthermore, 
as can be seen by the solid line in Fig. [3] 2(1bH c u/ A a is linear in l/d across the entire range of thicknesses. The 
dashed line in Fig. [3] represents the thickness dependence of the orbital critical field of Eq. Q using Eqs. JHJ and JJJ. 
Clearly, the critical field of bilayers with d < 5 nm is Zeeman mediated giving us a direct probe of the ISOC. 

In the context of Eq. one would attribute the linear critical field behavior in Fig.|2|to the thickness dependence 
coupling parameter b. In a related study Bergmann pil l25l | used magnetoresistance measurements to determine the 
SO scattering rate in non-superconducting Mg films coated with submonolayer coverages of Au. The Au produced 
large increases in an otherwise small intrinsic Mg SO scattering rate. In fact, the SO rate was observed to scale as 
b ~ l/d. This scaling behavior, however, is incompatible with the linear behavior in Fig.0 If one assumes that Eq. © 
is applicable for the thinnest bilayers, d < 5 nm, then the SO coupling parameter must scale as b ~ A /gP in order for 
the normalized critical field to be proportional to l/d. This suggests that the ISOC is not simply a scattering process. 
In fact, ISOC does not break time reversal symmetry, and is therefore not a pair-breaking interaction. Consequently, 
it is not meaningful to cast ISOC in terms os scattering processes, particularly in the T, H — > limit. Alternatively, 
it seems more appropriate to characterize ISOC as a boundary condition on the spin states of the superconductor. 
In a semi-infinite slab, for instance, the superconductor would have a non-negligible spin-triplet component near the 
Be-Au interface but decay into a pure spin singlet eigenstate as one moves into the interior. 

As discussed above, the spin susceptibility of the NCS bilayer is expected to be anisotropic in the strong ISOC 
limit. In order to address this issue we have measured the critical field of Be films and Be/Au bilayers as a function 
of tilt angle between the magnetic field and the plane of the film. Upon tilting the sample out of parallel orientation, 
orbital contributions to the critical field quickly dominate. Consequently, it is difficult to infer any anisotropy in the 
Zeeman response from the raw rotational data (see insets of Fig.|SJ). To circumvent this, we have measured the ratio of 
the critical field of Be/Au bilayers and Be films of equal thickness as a function of tilt angle 9. This, in effect, greatly 
improves the sensitivity to variations on the orbital background, thus giving us access to the angular dependence of 
the Zeeman contribution. Typical behavior is shown in Fig. 0] In order to better display the overall structure of 
the curves we have normalized the parallel field ratios to unity in this plot. Note the dip structure near 20° in each 
curve. We believe that this structure is consequence of an anisotropic susceptibility. In Fig. we plot the ratio of the 
critical field of a 5.4 nm Be/Au bilayer and a 5.4 nm Be film as a function of angle. The perpendicular critical field 
H C 2 of the bilayer was a factor of 3 higher than that of the Be film, hence the ratio is not unity at 6 = 90° . The solid 
line in Fig. [S]is the expected angular dependence of Eq. using an isotropic SO coupling parameter of b = 0.013 
mV for the Be film and b = 1.85 mV for the bilayer, as determined from fits to the inset data. For both samples 
A D ~ O.lmV. An isotropic b leads to a monotonic interpolation between the parallel and perpendicular field ratios. 
The isotropic solution is, in fact, monotonic independent of the relative thickness, T c , resistance, and/or perpendicular 
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critical fields of the two films. This reflects the fact that it is almost indistinguishable from what one would obtain 
using the phenomenological Tinkham formula 0,^3 which only involves H c ii and H C 2- In contrast, the dashed line 
depicts the solution to Eq. QJ assuming an exponentially attenuated SO parameter b = 1.85 exp(— 0/9 o ) mV, with a 
characteristic angle o = 2.5°. The dashed line captures the dip structure in the data which indicates that the Zeeman 
component of the critical field is angle dependent. The dip structure is consistent with an anisotropic superconducting 
susceptibility consisting of a large in-plane component (of the order of the normal state susceptibility) and a small 
transverse component. The sense of this anisotropy is of the opposite sign of that calculated for non-disordered 
2D NCS superconductors, where the in-plane component is expected to be half that of the transverse component. 
Nevertheless, the fact that ISOC induces a preferentially large in-plane susceptibility would also account for earlier 
reports of anomalously high parallel critical fields in Al/Pt bilayers [Tflj ]. 

In summary, we have used critical field measurements to probe the effects of interface spin-orbit coupling on the 
spin states of superconducting Be films. We find that a ~1 monolayer coating of Au on the Be surface produces a 
large, anisotropic enhancement to the Zeeman component of the critical field. The normalized parallel critical field 
H C \\/A Q of Be/Au bilayers scales as the inverse of the Be thickness and can exceed the Clogston limit by as much 
as an order of magnitude. The critical field scaling was observed over a wide range of A G which is inconsistent with 
the standard impurity scattering model of SO coupling. We believe that the scaling is, in fact, a manifestation of 
the superconductor's attempt to reconcile a mixed-spin boundary condition at the Au interface with the intrinsic 
spin-singlet ground state of Be. Naively, one would expect that the ISOC healing length would be of the order of 
£, but it may be significantly shorter in the presence of disorder. Nevertheless, the anisotropic spin susceptibility is 
clearly evident in the tilted field data. A high field study of the local tunneling density of states at the Be/Au and 
the Be/substrate interfaces, respectively, should provide an important local probe of the extent of spin- mixing at the 
two boundaries. 

We gratefully acknowledge enlightening discussions with Victor Edelstein, Gianluigi Catelani, Ilya Vekhter, Dana 
Browne, and David Young. This work was supported by the National Science Foundation under Grant DMR 02-04871. 
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FIG. 1: A) 0.1 x 0.1 /im atomic force micrograph of a 6 nm thick Be film evaporated onto glass at 84 K. B) Micrograph of 
Be film coated with 0.5 nm of Au. 
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FIG. 2: Transition temperature of Be/Au bilayers as a function of Be thickness. The Au thickness was 0.5 nm in each sampli 
The solid line is a best fit to the data using the empirical form of Eq. ©. Inset: relative diffusivity of Be/Au bilayers as 
function of Be thickness. The solid line is a fit to Eq. Q . 
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FIG. 3: Normalized parallel critical fields as a function of Be thickness for Be/Au bilayers (circles), Be/Pb (crossed box), pure 
Be films (triangles). The long dashed line represents the theoretical orbitally limited critical field given by Eq. Q. The solid 
line is a linear least fit to the bilayer data. The horizontal dashed line represents the Clogston critical field. 
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FIG. 4: Ratio of the critical field of Be/Au bilayers and Be films of equal Be thickness as a function of tilt angle, 9 = 
corresponds to parallel field, triangles: d = 7.4 nm, circles: d = 5.4 nm, diamonds: d = 2.9 nm. The solid lines are a guide to 
the eye. 
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FIG. 5: Ratio of the critical field of a 5.4 nm Be/Au bilayer (T c = 0.68 K, H c2 = 0.137 T, R = 240 Q) and 5.4 nm Be film 
(T c = 0.505 K, _f/ C 2 = 0.048 T, _R = 162 0) as a function of tilt angle. Solid line is the solution of Eq. (0 assuming an isotropic 
SO coupling parameter of b = 1.85 mV for the Be/Au bilayer and b = 0.0132 mV for the Be film. The dashed line is the 
solution of Eq. (0 assuming that the Be/Au SO parameter is exponentially attenuated with increasing tilt angle. Left inset: 
Fit of Eq. Q to the Be film data. Right inset: Solid line is a fit of Eq. Q to the Be/Au data with a constant b. The dashed 
line is a somewhat better fit using an exponentially form b = b exp(—9/0 o ). 



